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Analyzing Linear Communication Networks using 

the Ribosome Flow Model 

Yoram Zarai, Oz Mendel, and Michael Margaliot 


Abstract —The Ribosome Flow Model (RFM) describes the 
unidirectional movement of interacting particles along a one¬ 
dimensional chain of sites. As a site becomes fuller, the effective 
entry rate into this site decreases. The RFM has been used to 
model and analyze mRNA translation, a biological process in 
which ribosomes (the particles) move along the mRNA molecule 
(the chain), and decode the genetic information into proteins. 

Flere we propose the RFM as an analytical framework for 
modeling and analyzing linear communication networks. In this 
context, the moving particles are data-packets, the chain of sites 
is a one dimensional set of ordered buffers, and the decreasing 
entry rate to a fuller buffer represents a kind of decentralized 
backpressure flow control. For an RFM with homogeneous link 
capacities, we provide closed-form expressions for important net¬ 
work metrics including the throughput and end-to-end delay. We 
use these results to analyze the hop length and the transmission 
probability (in a contention access mode) that minimize the end- 
to-end delay in a multihop linear network, and provide closed- 
form expressions for the optimal parameter values. 

Index Terms —Packet-flow, communication networks, through¬ 
put maximization, end-to-end-delay, multihop, optimal hop- 
length, totally asymmetric simple exclusion process, statistical 
mechanics. 


I. Introduction 

The inflation of wireless devices and the dramatic increase 
in demand for real-time applications has led to a growing 
interest in the analysis and design of mobile ad hoc net¬ 
works (MANETs). Important features of MANETs include: 
nodes may operate as clients and as routers; flow control is 
decentralized; and the communication infrastructure between 
the nodes changes dynamically on the basis of network con¬ 
nectivity. The end-to-end features of such networks depend 
on the intricate correlations between the different flows in the 
system. The rigorous analysis of basic network characteristics, 
such as end-to-end delay, buffer utilization, and throughput, is 
a non-trivial challenge 

Wireless networks have been analyzed using different tools 
including information theory fTl and queuing theory ||2|, III- 
However, these often use simplified assumptions and yield ap¬ 
proximate results when characterizing the network throughput 
and delay performances. Eor example, large scale 2D networks 
have been studied based on Kleinrock’s Independence As¬ 
sumption a. According to this assumption, traffic correlations 
between different flows can be ignored and the effect on 
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delay performance is negligible, subject to sufficient traffic 
mixing and moderate-to-heavy traffic loads (see e.g., 0). 
However, recent studies suggest that this assumption does 
not hold in various real-world communication networks 0. 
Ad-hoc networks have also been modeled as a network of 
Jackson queues 0, Q. This is a stochastic model that 
assumes independent M/M/1 queues along the network. The 
service times are exponentially distributed and are mutually 
independent. As such, the average number of data-packets 
in each node is determined independently from the node’s 
neighbors and thus the interaction between adjacent nodes is 
neglected. 

Andrews et al. m point out that standard analysis tools 
may not be adequate for MANETs. They suggest several other 
tools, including non-equilibrium statistical mechanics. They 
state that: “In particular, the microscopic theories of vehicular 
traffic are promising since they do not treat vehicle flows as 
compressible fluids but explicitly focus on the dynamics of 
the individual vehicle. In other words, the theories treat the 
vehicular network as a system of interacting particles driven 
far from equilibrium, enabling the study of the dynamics of 
more general non-equilibrium systems-such as MANETs.” 

An interesting paper a (see also oni, m) has sug¬ 
gested using the Totally Asymmetric Simple Exclusion Pro¬ 
cess (TASEP) Ha to model and analyze ad hoc networks. 
TASEP is a stochastic model for particles that move along 
some kind of “tracks” or “trails”. A lattice of sites models 
the tracks. Each site may either be free or contain a single 
particle. Particles hop, with some probability, from one site 
to the consecutive one, but only if the target site is free 
(hence the term simple exclusion). This models particles that 
have volume and cannot overtake one another. In this way, 
the model captures the interaction between the particles. The 
term totally asymmetric is used to indicate that motion is 
unidirectional. 

TASEP and its variants are regarded as a paradigmatic 
model in non-equilibrium statistical mechanics, and have been 
used to model and analyze numerous natural and artificial 
multiagent systems, including traffic flow, molecular motors, 
mRNA translation and more KSi, m, in. Eor example, in 
the context of traffic flow, the particles [lattice] represent cars 
[the road]. 

Recently, Reuveni et al. ina studied a deterministic non¬ 
linear model, called the ribosome flow model (REM), that 
can also be derived by mean-field approximation of TASEP. 
Although nonlinear, the RFM is amenable to mathematical 
analysis using various tools including the analytical theory of 
continued fractions 1(161 . ifTTl . monotone systems theory ifTSl . 
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lfT9l . contraction theory Il20l . convex optimization theory 11211 . 
and matrix theory ll22l . 

We propose the RFM as a theoretical framework for linear 
communication networks. In this context, the moving particles 
are data-packets and the chain of sites is a one dimensional 
set of ordered buffers. For the RFM with homogeneous link 
capacities, we provide closed-form expressions for important 
features including: steady-state buffer occupancy, network 
throughput, and end-to-end delay. We use these results to 
provide closed-form expressions for the hop length and the 
transmission probability that minimize the end-to-end delay 
in ad-hoc multihop linear networks. 

The remainder of this paper is organized as follows. Sec¬ 
tion [n] reviews the RFM. Section |III] explains how the RFM 
can be used to model a linear communication network. Sec¬ 
tion presents the main results. Section 0 describes two 
applications of these results to a multihop network. Section IVll 
concludes and describes possible directions for future research. 
To streamline the presentation, all the proofs are placed in the 
Appendix. 

II. Ribosome Flow Model 

The RFM is a set of n nonlinear first-order ordinary 
differential equations (ODEs): 

±1 = Ao(l - xi) - - X2) 

±2 = XlXl{l - X 2 ) - X 2 X 2 {^ - X 3 ) (1) 

Xn — An—la^n—1(1 Xn') ^nXn- 

This is a compartmental model ll23l . consisting of n compart¬ 
ments (or sites) ordered along a linear ID chain. Each site 
may contain some material. (In the case studied here, each 
site is a buffer that contains data-packets of fixed length.) The 
state-variable Xi{f), i G {1,... ,n}, represents the occupancy 
level of site i at time t. The occupancy levels are normalized 
so that Xi{t) G [0,1] for all t > 0, with Xi{t) = 0 [Xi(t) = 1] 
representing that site i is completely empty [full] at time t. 
Site 1 is fed by a source with entry rate Aq > 0. Material flow 
from site i to site i-\-l depends on a elongation rate Xi > 0. 
The material exits at the end of the chain, that is, from buffer n 
with exit rate A„ > 0. 

To explain the equations of this model, consider the equation 
for xi, i.e., the rate of change in occupancy level at site 1. 
The term Ao(l — a;i) is the transition rate into the chain: as Xi 
increases this rate decreases and, in particular, when xi = 
1 (i.e., buffer 1 is completely full) this rate becomes zero. 
The term Aia;i(l — X2) is the transition rate of material from 
site 1 to the consecutive site 2. This is proportional to the 
occupancy level at site 1, and also to (1 — X 2 )- In particular, 
when a :2 = 1 this becomes zero. Thus, as site 2 becomes 
fuller the flow from site 1 to site 2 decreases. In this way, 
the REM encapsulates both the simple exclusion principle and 
the unidirectional movement of TASER The exit rate from the 
last site is R{f) := XnXnit). (In the biological context, this 
represents the protein production rate.) 


The state space of the REM is the n-dimensional unit 
cube C" := [0,1]". Let int(C"') denote the interior of C”, 
and let x{t,a) denote the solution of the REM at time t for 
the initial condition x{Q) = a. 

Proposition 1 MSy The RFM admits a unique equilibrium 
point e = e(Ao,...,A„) G int(C”), and for any a G C", 
x(t, a) G int((7") for all t > 0, and limi_>.oo x(t, tco) = e. 

This means that every trajectory of the REM converges to a 
steady-state that depends on the rates Aq, ..., A„. In particular, 
R{f) converges to the steady-state exit rate R := Xn^n- 

The simulation results reported in lITSll show that REM 
and TASEP yield similar predictions of steady-state rates. 
However, recent research results on the REM suggest that 
expressions for important quantities are considerably simpler 
than the corresponding expressions in TASEP. 

Several recent papers analyzed the REM. In the particular 
case where Ai = • • • = A„ := Ac, with Ac denoting 
the common value, the REM is called the Homogeneous 
Ribosome Elow Model (HREM). This model includes only 
two parameters Aq and Ac, so its analysis becomes more 
tractable ca. The HREM with an infinitely-long chain, (i.e. 
with n 00 ) was considered in iflTll . There, a simple closed- 
form expression for Coo := limn^oo e„ was derived, as well as 
explicit bounds for |eoo—e„ | for all n > 2. A closed-loop REM 
(representing mRNA circularization) has been studied in |[T9l . 
It has been shown that the closed-loop system admits a unique 
globally asymptotically stable equilibrium point. Ref. Il20l has 
shown that when the rates Xi are periodic functions of time, 
with a common period T, then ([T]i admits a unique periodic 
solution p(t), with period T, and x(t,xo) — p{t) for every 
initial condition xq G C". In other words, the REM entrains 
to periodic excitations in the rates. Ref. im has shown that 
the steady-state exit rate R = i?(Ao, •.., A„) is a strictly 
concave function of Ao,...,A„ on This means that 

the problem of maximizing the steady-state exit rate, under an 
affine constraint on the rates, is a convex optimization problem. 

HI. Modeling Linear Communication Networks 
USING THE REM 

Fig. □ illustrates the general block diagram of a linear 
communication network with multiple relay nodes modeled 
using the REM. Here the material flowing is data packets, 
each site corresponds to a communication node that contains 
a buffer, and Xi (f) is the (normalized) occupancy of buffer i 
at time t. The parameter Ai > 0 is the capacity of the link 
from node i to node i-\-l. The network throughput at time t 
is R[t). The equation Xi = Xi-iXi-i{l — Xi) — XiXi{l — Xi+i) 
represents a decentralized flow control algorithm, as the flow 
through buffer i depends on the occupancy at buffers i — 1 ,i, 
and i-\-l only. 

Prop. [U means that for any set of strictly positive link 
capacities the flow control algorithm is stable in a strong 
sense, as the dynamics keeps every buffer from overflow or 
underflow at all time t (QoS flows are one example where 
buffer overflow must be avoided). Eurthermore, the network 
always converges to a steady state in which the occupancy of 
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buffer i is ei G (0,1). Note that may also be interpreted as 
the probability that site i is occupied at steady-state. 

Recall that the backpressure algorithm (see, e.g., 1251, 

f26i . Il27l . schedules packets from node i to node i + 1 
based on the differential backlog max{0,Xi — Xi+i}. As Xi 
increases and Xi+i decreases, the effective flow rate increases. 
This also happens in the RFM, as the flow rate is XiXi{l — 
Xi+l). 

Note that TASEP is a stochastic, discrete-time, asyn¬ 
chronous model, whereas the RFM is a deterministic, 
continuous-time, synchronous model. Indeed, in the RFM all 
the XiS are updated synchronously. Conceptually, the RFM 
[TASEP] is more suitable for modeling access protocols like 
slotted-ALOHA or CDMA [TDMA] where the nodes [do not] 
transmit at the same time. 

In the context of communication networks, TASEP corre¬ 
sponds to a chain of buffers, where each buffer can either 
contain a packet or not. The RFM provides a more flexible 
framework, as the buffer level Xi takes values in [ 0 , 1 ], so 
the buffer can be either empty (xi = 0 ), full {xi = 1 ), half 
full (xi = 1/2), and so on. As shown in Section HVl analysis 
of the RFM leads to simpler expressions than those known 
for TASEP. Nevertheless, below we also simulate the time- 
averaged steady-state behavior of TASEP and show that the 
results agree well with the analytical formulas derived for 
the RFM. 


IV. Analysis of Linear Communication Networks 
USING THE RFM 

In this section, we analyze properties of the RFM that are 
relevant for communication networks. We consider the special 
case where all the link capacities are equal, i.e. Aq = • • • = 
An := Ac. We refer to this case as the totally homogeneous 
RFM (THRFM). 


Proposition 2 Consider the THRFM with n sites. The steady- 
state occupancy level at site i is 


1 1 
ei = 2 + 2 


cot 


-f 3 

The steady-state throughput is 
R = 


{i -b l)7r 

n -b 3 


4cos2(;j^)' 


i = 1,..., n. 

( 2 ) 

(3) 


The steady-state delay experienced by a data-packet at node i 
is 

2 


Di = — cos 
Ar 


n -b 3 


and the end-to-end delay is 


De2e ■=J2^^ = T- 


, TT , ,(i + 1 )t A 

TL o Ti o J 

(4) 

cos" I -r I ■ (5) 


2 n 


n -b 3 


Eq. (|2]i implies that the steady-state buffer occupancies are 
independent of the homogeneous link capacity Ac (but Ac 
does control the convergence rate to the steady-state). Com¬ 
bining dU with the identity cot (a:) = — cot(7r — a:) implies the 
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Fig. 2. Analytical (O) and empirical (+) steady-state occupancy level as 
a function of 2 = 1,, 19, in a THRFM with n = 19 nodes. 


symmetry relation et ~ | and this gives 

n 

'^ei=^nl2. ( 6 ) 

This implies that for n odd, e(n+i )/2 = 1/2. Eq. ^ also im¬ 
plies that ei > • • • > e„. This means that at steady-state, much 
of the queuing occurs closer to the source node, and preference 
is given to nodes closer to the destination node. Indeed, it 
was suggested that a scheduling algorithm giving preference 
to serving links closer to the destination node achieves optimal 
end-to-end buffer usage in linear networks 1291 . 

Eig. I 2 ] depicts the e^s for a network with n = 19 nodes. It 
may be seen that decreases with i. Note also the symmetry 
around eio = 1/2. Eor example, eg — 5 = | — 612 - Eig. |2]also 
depicts empirical occupancy levels obtained via simulation of 
TASEP with n = 19 nodesQ 

Eq. m implies that: (1) the steady-state throughput R is 
proportional to the homogeneous link capacity Ac; (2) R is 
upper-bounded by Ac/2 for all n; and (3) R decreases with 
the network length n, with lim„_>oo R = Ac/4 (see Eig. [3). 

Eq. ® implies that Di > ■ ■ ■ > Dn, as the steady-state 
delay experienced by a data-packet at node i is proportional 
to the occupancy level at that node. The end-to-end delay De 2 e 
is inversely proportional to the link capacity Ac, and increases 
more or less linearly with the network length n. In partic¬ 
ular, De 2 e ~ 2nl\c for large values of n. Eig. |4] depicts 
De 2 e as a function of Ac for a THREM with n = 19. It also 
shows the empirical values of De. 2 e obtained via simulating 
TASEP and averaging the end-to-end delays of data-packets 
(after discarding the first 350, 000 cycles). 


*The TASEP simulation used a parallel update mode. At each time tick tj., 
the sites are scanned from site n backwards to site 1. If it is time to hop and 
the consecutive site is empty then the particle advances. If the consecutive 
site is occupied the next hopping time, ij. -l-ej., is generated. Eor site i, is 
exponentially distributed with parameter (1/Ai). The occupancy at each site 
is averaged through the simulation, with the first 350, 000 cycles discarded 
in order to obtain the steady-state value. 
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Fig. 1. Linear communication network with n relay-nodes modeled using the RFM. SRC [DST] represents the source [destination] node. The circles represent 
the relays. State-variable Xi{t) G [0,1] is the normalized number of data-packets at buffer i at time t. The capacity of the link from site i to site 2 + 1 is Aj. 
The throughput at time t is R(t) = XnXnit). 
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Fig. 3. Throughput R for link capacity Ac = 1 as a function of n. 



Fig. 4. Steady-state end-to-end delay De 2 e as a function of the link 
capacity Ac in a network with n = 19 nodes. 


V. Application to Multihop Networks 

Ad-hoc multihop networks and their military and civilian 
applications are attracting considerable interest, yet there 
remain significant theoretical challenges in the performance 
analysis of these networks 1301, El, il, and thus in their sys¬ 
tematic design. The decentralized, multiagent network struc¬ 
ture and the possibility for multihop routing make the analysis 
of fundamental communication metrics, including throughput, 
latency and capacity quite challenging. In multihop networks. 


data-packets may be routed over many short hops or over a 
smaller number of long hops. An important research topic is 
determining various parameters, e.g. the hop length, that opti¬ 
mize some performance metric, e.g. minimize the end-to-end 
delay. (Of course, in general there could be additional factors 
impacting the decision between short and long hops HD.) 
In this section we show how to use the analysis results on 
the THRFM described above to systematically address this 
type of optimization problems. 

We assume that all the nodes in the network use the same 
channel. The nodes are equally spaced along the communi¬ 
cation system, with a separation £ between any two adjacent 
nodes. Attenuation in the channel is modeled as the product 
of a Rayleigh fading component and a large-scale path loss 
component with exponent 7. The channel noise is taken 
to be additive white Gaussian noise (AWGN) with power 
spectral density Nq. We define the transmission from node i 
to node f -|- 1 to be successful if the (instantaneous) signal-to- 
interference-and-noise ratio (SINR) at node i-l-1 is greater than 
a predetermined threshold 0. The probability of a successful 
transmission is thus ps := Pr(SINR > 0). 

1) Optimal Hop Length at SNR Limit Regime: Consider an 
ad-hoc multihop linear network where the main effect on a 
successful transmission is the channel attenuation and noise. 
This happens, for example, when the nodes are relatively far 
so that the effect of the spatial reuse is negligible relative to 
the channel noise. In this regime, 

Ps = Pr(SNR > 0) = exp(-0Vor). (7) 

We model this using the THRFM with link capacities Ac = 
rps, where r > 0. For simplicity, we take r = 1. The end-to- 
end delay in the THRFM given in (|5]) then becomes 

£>e2e = 2nexp(0jVor) cos^ ( ^ 

\n + 6j 

The closed-form expression for De 2 e allows us to address 
various optimization problems. For example, consider the 
problem of finding the hop length that minimizes the end- 
to-end delay. Let ml denote a hop of m steps of unit length 1. 
Replacing I with ml and n with n/m in ® yields 

De 2 e = — exp(0Vo(mO^) COS^ ( , / . q ) ■ (9) 

m \(n/m) + 3/ 
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Fig. 5. Optimal hop length m* as a function of the SNR threshold © and 
the path loss exponent 7 for a network with n = 200 nodes. 


Fig. 6 . Optimal transmission probability q* as a function of the threshold © 
and the path loss exponent 7 . 


In particular, for n/m ^ 1, 

— Stt- 

De2e ~ De2e '■= - exp(0iVo (mZ)^) . (10) 

m 

This expression shows that m affects Z?e 2 e in two different 
ways. On the one-hand, increasing m decreases the total 
length of the chain and thus decreases the delay. On the 
other-hand, increasing m increases the separation between 
the communicating nodes. This decreases the probability of 
a successful transmission thus increasing the delay. 

Proposition 3 The hop length that minimizes the end-to-end 
delay De 2 e in the SNR limit regime is either \m*~\ or 
where 

If 1 

and for this value, the optimal end-to-end delay is 

D*e 2 e = ‘^nl ( 70 iVo)^/'^ exp(l/ 7 ). 

Eq. (E) implies that m* is monotonically decreasing with 
0 and Nq. This is reasonable, as an increase in any of 
these parameters decreases the effective SNR at the receiving 
node, which in turn decreases the probability of a successful 
transmission. This can be compensated by decreasing the hop 
length. Fig. [^depicts the analytical and empirical values of m* 
as a function of 0 and 7 for f = 1 and Nq = —50 dBW/Hz. 
The empirical values were obtained by simulating several 
TASEPs, each one with 200/m sites and with an exponentially 
distributed next time to hop with parameter Ac = Ps as given 
in ©, but with i replaced by mi. 

2) Optimal Transmission Probability at SIR Limit Regime: 
Consider now an ad-hoc multihop wireless network where the 
main effect on a successful transmission is the transmission 
interference. Let q denote the contention probability, i.e. the 
probability that a node transmits at each time slot. In this case, 

Ps exp(- 9 c/ 2 ), (12) 

X 

where c := “ 1 (^^e 1331). 


A natural question is how to determine the contention 
probability that minimizes De 2 e- To address this, we model 
this scenario as a THRFM with Ac = qps- Substituting this 
in 0 yields 

De 2 e = {v/q)exp{qc/2), (13) 

where v := 2ncos^ This expression shows that q 

affects De 2 e in two different ways. On the one-hand, increas¬ 
ing q increases the link capacity and thus decreases the delay. 
On the other-hand, increasing q increases the interference. This 
decreases the probability of a successful transmission thus 
increasing the delay. 

Proposition 4 The contention probability that minimizes the 
end-to-end delay in the SIR limit regime is 

q* = min{l, 2/c}. (14) 

Fig. | 6 ] depicts the analytical and empirical optimal con¬ 
tention probability q* as a function of 0 and 7 . Note that q* 
is inversely proportional to both 0 and to 7 . 

VI. Discussion 

TASEP is an important model in non-equilibrium statistical 
mechanics that has been used to study numerous biological and 
artificial multiagent systems including multihop wireless net¬ 
works El, m, im. The REM is a “deterministic substitute” 
of TASEP that is highly amenable to mathematical analysis. 
In this paper, we proposed the REM as a new theoretical 
framework for modeling and analyzing linear communication 
networks. This has several advantages. The dynamics of the 
REM guarantees that the buffer levels do not overflow or 
underflow for all time t, and converge to a unique steady-state 
level. For the particular case of the THRFM, another important 
advantage is that simple closed-form expressions for important 
network metrics such as the steady-state throughput and end- 
to-end-delay exist. These expressions can be used to select 
optimal parameters, for example, the optimal hop length in a 
multihop network. 
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For the more general case of an RFM with non- 
homogeneous rates closed-form expressions for the network 
metrics are not known. However, they can be easily computed 
in an efficient and numerically stable way. Indeed, it has been 


shown 

in EB 

that the 

(n -f 2 ) 

X (n 

- 1 - 2 ) matrix 
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,- 1/2 
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-1/2 
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0 A 

- 1/2 

n 
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(15) 

has real and distinct eigenvalues: Ci 

< C 2 < • • • 

< Cn-l-2 

and 


that the steady-state throughput in the RFM is i? = 

This means that R, and then every other steady-state metric, 
can be determined using numerical algorithms for calculating 
the eigenvalues of symmetric (tridiagonal) matrices. 

The link capacities in a communication network must often 
be constrained, e.g. to reduce interference or because of energy 
constraints 041 . A natural question is how to maximize the 
throughput subject to a constraint on the capacities. In the 
linear communication network modeled using the RFM, this 
leads to the following optimization problem; 

Problem 1 Fix b,wo,wi,... ,Wn > 0. Maximize R = 
i?(Ao,..., A„), with respect to its parameters Ao,...,Ati, 
subject to the constraints: 

n 

< b, ( 16 ) 

Aq; ■ • ■ 5 ^ 0- 

The values Wi can be used to provide different weighting to the 
different link capacities. It has been shown in OTI that i? is a 
strictly concave function of the A^s, and since the constraint on 
the rates is affine, ProbIem[T]is a convex optimization problem. 
Thus, the solution A*, i = 0,..., n, is unique and can be found 
using efficient numerical algorithms that scale well with n. 

Fig. |7] depicts the optimal link capacities A*, j = 0,..., 9, 
for a communication network with 9 relay nodes and the 
constraint (fTSl) with wq = • • • = Wn = 6 = 1 . It may be 
seen that the A*s are symmetric, i.e. A* = Ag_j, and that they 
increase towards the center of the network. In other words, the 
capacities near the center of the network have a higher effect 
on the throughput than those near the ends of the network. 

Another interesting research topic is to model and analyze 
more complex network topologies, comprising of multiple 
flows that share common nodes, using a set of interconnected 
RFMs. 
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Fig. 7. Optimal link capacities A* as a function of i for a linear 
communication network with n = 9 nodes. 

Appendix: Proofs 

Proof of Prop. H It has been shown in ISl that for 
an HRFM with m sites and Aq —oo the equilibrium 
point e = [ei,...,em] satisfies e* = n ^ ) , i = 

l,...,m, where d^ := cos(;;^)-f sin(^) cot(^). In 
particular, ei = 1. Since Ai = • • • = A„ = Ac, ([B 
implies that the steady-state satisfies [ei 62 • • • Cm] = 

[1 ei • • • Cm-i], where [ei • • • Cm-i] is the equilib¬ 
rium of a THRFM with dimension m — 1. Taking m = n + 1 
yields (|2]l. At steady-state, the exit rate is i? = XnCn = XcCn, 
and this is equal to the rate of packet flow across each node. 
Combining this with ^ yields ([B- By Little’s Theorem ||3, 
the delay at node i at steady-state is Di = Ci/R. Combining 
this with 0 and 0 yields (|4|i. The end-to-end delay is 
ELi = (1/-R) ELi using ® and 0 yields 0 . ■ 

Proof of Prop. [3] Differentiating Z 3 e 2 e with respect to m 
yields ^^Z)e 2 e = ^ exp{Noe{lmr ) {-/NoQ^lmr - 1 ), 
= |^exp(z)( 72 z 2 + ^:^(^_3) + 2 ), where z := 
. It is straightforward to verify that -^^De 2 e > 0 
for all m > 0 and 7 > 1 , thus for all practical cases 
De 2 e is a strictly convex function of m, and the unique 
global minimum m* can be found by equating -^De 2 e to 
zero, which yields (HB (of course, the relevant value is 
either \m*~\ or as the hop length must be an integer). 

Substituting m* in (fTol i completes the proof. ■ 

Proof of Prop. [3 Differentiating (fT^ yields -^De 2 e = 
z;exp(cg/ 2 )^^, and ■§^De 2 e = wexp(c(?/ 2 )^^, where 
f{z) := — 4:Z + 8 . Since v > 0 and f(z) > 0 for all z, it 

follows that -^De 2 e > 0 for all g > 0 , so De 2 e IS a strictly 
convex function of q, and the unique global minimum q* can 
be found by equating ^fAe 2 e to zero. This yields q* = ^jc. 
If c > 2 then 5 * = ^jc < 1, so this is a valid minimizer (recall 
that q S (0, 1]). If c < 2 then 2lc> 1, so the valid minimum 
is q* = l.m 
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